Origin of the X-ray emission in the nuclei of FRIs 
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Abstract. We investigate the X-ray origin in FRIs using the multi- waveband high resolution 
data of eight FR I sources, which have very low Eddington ratios. We fit their multi-waveband 
spectrum using a coupled accretion-jet model. We find that X-ray emission in the source with 
the highest Lx (~ 1-8 x 10~''LEdd) is from the advection-dominated accretion fiow (ADAF). 
Four sources with moderate Lx (~ several xlO~^I/Edd) are complicated. The X-ray emission of 
one FRI is from the jet, and the other three is from the sum of the jet and ADAF. The X-ray 
emission in the three least luminous sources (Lx < 1.0 x 10~^LEdd) is dominated by the jet. 
These results roughly support the predictions of Yuan and Cui (2005) where they predict that 
when the X-ray luminosity of the system is below a critical value, the X-radiation will not be 
dominated by the emission from the ADAF any longer, but by the jet. We also find that the 
accretion rates in four sources must be higher than the Bondi rates, which implies that other 
fuel supply (e.g., stellar winds) inside the Bondi radius should be important. 



1. Introduction 

Radio galaxies are usually classified as FRI or FRII sources depending on their radio 
morphology. FRI radio galaxies (defined by edge-darkened radio structure) have lower radio 
power than FRII galaxies (defined by edge-brightened radio structure due to compact jet 
terminating hot spots) (FanarofF & Riley 1974). Different explanations of division of FRI and 
FRII radio galaxies invoke either the interaction of the jet with the ambient medium or the 
intrinsic nuclei properties of accretion and jet formation processes (e.g., Bicknell 1995; Reynolds 
et al 1996a; Hardcastle et al 2007). Accretion mode in low power FRIs may be different from 
that in powerful FRIIs. There is growing evidence to suggest that FRIs type radio galaxy nuclei 
possess advection-dominated accretion flow (ADAF; or "radiative inefficient accretion flows"; see 
Narayan & McMlintock 2008 for a review). In fact, we now have strong observational evidence 
that ADAFs may be powering various types of low-luminosity active galactic nuclei (LLAGNs; 
e.g., Yuan 2007; Ho 2008 for reviews), not only FRIs. It was found that the hard X-ray photon 
indices of both X-ray binaries (XRBs) and AGNs are anti-correlated with the Eddington ratios 
when the Eddington ratio is less than a critical value, while they become positively correlated 
when the Eddington ratios higher than the critical value (e.g., Wu and Gu 2008). This results 
provide evidence for the accretion mode transition near the critical Eddington ratio. 

One of the uncertainties in FR Is is the respective contribution of ADAFs and jets at various 
wavebands. The least controversial nuclear emission are the radio and also optical emission, 
which is believed to be dominated by the jet (e.g., Wu and Cao 2005; Chiaberge et al 2006). It 
is still an open question for the origin of the core X-ray emission in FRIs (and also LLAGNs). 



One possibility for the X-ray emission is dominated by the ADAF (e.g., Reynolds et al. 1996b; 
Quataert et al 1999; Merloni 2003). The other possibility is come from the jet (e.g., Falcke et al 
2004; Markoff et al 2008), and recent observation on some individual extremely low luminosity 
sources supported this possibility (e.g., Garcia et al. 2005 for M31; Pellegrini et al 2007 for 
NGC821, etc.). Then an important question is systematically in what kind of condition the 
radiation from the jet or ADAF will be important in X-ray band. 

2. Sample 

The FRI sample used for the present investigation is selected from Donate et al (2004), which 
have the estimated black hole mass, Bondi accretion rate, optical, radio, and X-ray nuclear 
emission (Tabic 1). There are 9 FRIs in their sample which have compact core X-ray emission 
and have been observed by Chandra. We excluded 3C 270 since the optical emission may be 
obscured by its large intrinsic column density {Nh ~ lO^^cm"^). Therefore, our final sample 
include 8 FRIs. 

3. Coupled accretion-jet model 

The accretion flow is described by a hot, optically thin, geometrically thick advection-dominated 
accretion flow. Following the proposal of Blandford and Begelman (1999), we parameterize the 
radial variation of the accretion rate with the parameter p„ caused by the possible wind (e.g.. 
Stone et al 1999), M = Mout(-R/-Rout)^") where Mout is the accretion rate at the outer boundary 
of the ADAF -Rout- We calculate the global solution of the ADAF. The viscosity parameter a and 
magnetic parameter P (defined as ratio of gas to total pressure in the accretion flow, (3 = Pg/Ptot) 
are fixed to be a = 0.3 and f3 = 0.9. Another parameter is 6, describing the fraction of the 
turbulent dissipation which directly heats electrons. Following Yuan et al (2006), we use S = 0.3 
and = 0.25 in all our calculations. The radiative processes we consider include synchrotron, 
bremsstrahlung and their Comptonization. We set the outer boundary of the ADAF at the 

Bondi radius Rb = 2GMbh/c5, where Cg = jkT / iinip is the adiabatic sound speed of the gas 
at the Bondi accretion radius, T is the gas temperature at that radius, f^ = 0.62 is the mean 
atomic weight, rup is the proton mass, and 7 = 4/3 is adiabatic index of the X-ray emitting gas. 
After the ADAF structure is obtained, the spectrum of the flow can be calculated (Yuan et al 
2003). 

The jet model adopted in the present paper is based on the internal shock scenario, widely 
used in interpreting gamma-ray burst (ORB) afterglows (e.g., Piran 1999). A fraction of the 
material in the accretion flow is assumed to be transferred into the vertical direction to form a 
jet due to the possible standing shock near the black hole (BH). From the shock jump conditions, 
we calculate the properties of the post-shock flow, such as electron temperature Tg. The jet is 
assumed to have a conical geometry with half-open angle (p and bulk Lorentz factor Tj which 
are independent of the distance from the BH. The internal shock in the jet should occur as a 
result of the collision of shells with different Tj , and these shocks accelerate a small fraction of 
the electrons into a power-law energy distribution with index p. We assume that the fraction 
of accelerated electrons in the shock = 10% and a typical value of = 0.1 (e.g., Laing & 
Bridle 2002) in our calculations. Following the widely adopted approach in the study of GRBs, 
the energy density of accelerated electrons and amplified magnetic field are described by two 
free parameters, ec and ee- Obviously, and ec are not independent. The Lorentz factor of 
most FRIs are determined from the observations and a typical value Tj = 2.3 (corresponding 
to v/c = 0.9, e.g., Verdoes Kleijn et al 2002; Laing &; Bridle 2002) is adopted for those not 
having estimations, which will not affect our main conclusion since that our results are not very 
sensitive to the Lorentz factor. 



4. Results 

We use the ADAF-jet model to fit the spectrum of FRIs. As we state in Sect. 1, the radio 
emission, and perhaps optical as well, is from the jet. Based on the assumption to the jet model 
described in Sect. 3, the contribution of the jet to the X-ray band is well constrained once we 
require the jet model to fit the radio and optical data. We then adjust the parameter of the 
ADAF and combine it with the jet contribution to fit the X-ray spectrum. 

4.1. ADAF dominated the nuclear X-ray emission for higher Eddington ratio sources 
The radio morphology and power of 3C 346 would rank as either a low-power FR II source or a 
high-power FRI. Chandra observations have detected an unresolved core with Lx(2 — 10 keV) = 
1.9 X lO''^ erg s'^ and a photon index of F = 1.69±0.09 (Donato et al 2004). Figure 1 shows the 
fitting result. The parameters of the jet are rhjet = 3.5 x 10~^, eg = 0.14, €b = 0.02, and p = 2.4. 
We find that the jet model can describe well the radio and optical data. But the X-ray emission 
of the jet model is several times lower than the Chandra observations, and the X-ray data can 
be well fitted by the ADAF. The required accretion rate is rhout = 2.8 x 10~^. The ratio of mass 
loss rate in the jet to accretion rate of ADAF at 10 Rs is about mjet/"^(10-Rs) =0.9%, where 
i?S is the Schwarzschild radius. This source is relatively luminous, with Lx/i^Edd = 1-8 x 10"^, 
and the X-ray emission is dominated by the ADAF, not by the jet. 




Figure 1. Spectral energy distribution Figure 2. The same as Fig. 1, but for 

of 3C346. The dot-dashed, dashed, and 3C31. 

the solid lines show the emissions from the 

ADAF, jet, and their sum, respectively. 

The thin long-dashed line is synchrotron- 

self-Compton spectrum of the jet. 



4.2. ADAF+jet for moderate Eddington ratio sources 

3C 31 is a twin-jet FRI radio galaxy. The X-ray spectrum of the core is quite flat, with a photon 
index F = lAStoH and Lx(2 - 10 keV) = 4.7 x lO^'Ocrg s"! (~ 4.4 x lO-^Ledd) (Evans et al 
2006). Figure 2 shows the fitting result of 3C 31. The parameters of the jet are mjet = 2.7x 10~^, 
ee = 0.2, ee = 0.02, and p = 2.5. We find that the radio, optical and even the soft X-rays nuclear 
emission (e.g., 1 keV) can be well-fitted by a pure jet model. However, the hard X-rays cannot be 



fitted by a jet, but can be well-fitted by the ADAF with an accretion rate of rhout = 3.7 x 10~^. 
The ratio, rhjet/fniW Rs), is about 9%. We can see that jet contribution is important at soft 
X-ray band, while the ADAF contribution is important at the hard X-ray band. 

This is also the case for other FR Is with moderate Eddington ratios {Lx / L^dd ~ several x 
10"^, e.g., 3C317, B2 0055-^30, and possibly also 3C449), except B2 0755-^37, where all the 
emission can be fitted by a pure jet model (see Wu et al 2007 for more details). 

4-3. Jet dominated the nuclear X-ray emission for lower Eddington ratio sources 
The point-like nucleus is detected in 3C 66B by Chandra yielding Lx(2 — 10 keV) = 1.1 x 
lO^^ergs"^ (~ 1 X lO-^LEdd) and a photon index of F = 2.171^;^^ (e.g., Donato et al 2004). 
Figure 3 shows the fitting result of nucleus of 3C 66B. We find that all the radio, sub-millimetre, 
optical and X-ray can be fitted by a pure jet model very well (dashed line). The parameters 
of the jet are mjet = 1 x 10~^, eg = 0.18, cb = 0.02, and p = 2.35. For illustration purposes, 
we also show the X-ray emission using an ADAF model with rriout = 2.6 x 10~^ (dot-dashed 
line). The predicted spectrum by an ADAF is too hard to be consistent with the observations 
(so the accretion rate in the ADAF should be smaller than 2.6 x 10~^). The X-ray emission in 
this source should be produced by the jet. 

3C 272.1 is the lowest Eddington source (~ 6.8 x 10~^LEdd)) for which many waveband 
high resolution observations from radio to X-ray have been carried out. Figure 4 shows the 
fitting result of nucleus of 3C 272.1. The parameters of the jet are mjet = 4.9 x 10~^, eg = 0.28, 
ee = 0.005, andp = 2.5. We can sec that the radio, sub-millimetre, optical, and especially, the X- 
ray emission can be fitted by the jet model very well. On the other hand, the predicted spectrum 
by an ADAF (with rhout = 1-9 x 10~^) is too hard to be consistent with the observations. So in 
this source, the X-ray emission is also dominated by the jet. 
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Figure 3. The same as Fig. 1, but for 3C Figure 4. The same as Fig. 1, but for 3C 
66B. 272.1. 



4.4- Fuel supply in these FRIs 

Typically, the Bondi accretion rate is a good estimation to the mass accretion rate. However, 
Pellegrini (2005) show that there is no relation between the nuclear X-ray luminosity and the 
Bondi accretion rate in LLAGNs, and the X-ray emission of some sources is higher than the 



values predicted by ADAFs with the Bondi accretion rate. The Bondi accretion rate in our 
sample has been estimated (see, Donato et al 2004). We find that the accretion rates mout 
required in our model of four FRIs (3C 346, 3C 31, 3C 449, and 3C 317), among the five in 
which we can have good constraints to their accretion rates, are higher than their Bondi rates 
by factors of 9, 18, 112, and 1.05, respectively (see also table 1). Given that the radial velocity 
of the accretion flow is acs , a more accurate estimation of the accretion rate is mout ~ CK?7^Bondi 
where a is the viscous parameter (Narayan 2002). Therefore, the Bondi accretion rate is only 
a lower limit of the real rate and other fuel supply must be important, such as the gas released 
by the stellar population inside the Bondi radius (Soria et al 2006; Pellegrini 2007). 



Table 1. Accretion and jet properties. 



Smircp Inp-inCMBHN Lx(2-10keV) rhj^t mCRo'l" ^Kin 

source ioglU( ) ^^^^ rh(iORs) ^K^b) rn^ m{ioRs)c^ 



3C 346 


8.89 


1.8 


X 


10- 


-4 


0.91% 


2.8 X 10^2 


3.1 


X 


10- 


-3 


0.03 


B2 0755+37 


8.93 


5.2 


X 


10" 


-6 


> 3.9% 


< 8.9 X 10-3 


3.2 


X 


10" 


-2 


> 0.10 


3C 31 


7.89 


4.4 


X 


10" 


-6 


9.0% 


3.7 X 10-3 


2.1 


X 


10^ 


-4 


0.19 


3C 317 


8.80 


3.4 


X 


10" 


-6 


8.9% 


4.7 X 10-3 


4.5 


X 


10" 


-3 


0.28 


B2 0055+30 


9.18 


2.4 


X 


10" 


-6 


3.5% 


2.7 X 10-3 


1.4 


X 


10" 


-2 


0.08 


3C 66B 


8.84 


1.0 


X 


10- 


-6 


> 5.9% 


< 2.6 X 10-3 


2.5 


X 


10- 


-2 


> 0.18 


3C 449 


8.42 


8.0 


X 


10- 


-7 


14.3% 


1.9 X 10-3 


1.7 


X 


10- 


-5 


0.44 


3C 272.1 


8.35 


8.3 


X 


10- 


-8 


> 7.3% 


< 1.9 X 10-3 


6.0 


X 


10" 


-2 


> 0.22 



(a) m(i?B) is the dimensionless accretion rate at the Bondi radius through our spectra fitting; 

(b) mB is the dimensionless Bondi accretion rate estimated from the X-ray observation. 



5. Discussion and Conclusion 

Two possibilities exist for the X-ray origin in FRIs: ADAF or jet. We use a coupled ADAF-jet 

model to fit the multi-waveband spectrum to try to investigate this problem. We find that 
the jet can well describe the radio and optical spectra for most FRIs in our sample. The soft 
X-ray flux at ~lkeV of all FRIs is roughly consistent with the predictions of the jet. This 
result indicates why a tight correlation is found among radio, optical and soft X-ray (Evans et 
al 2006; Balmavcrdc et al 2006). For 3C 346, the highest X-ray luminosity source in our sample 
which has Lx = 1.8 x lO-^Lgdcl) its X-ray spectrum is dominated by the ADAF and the jet 
contribution is negligible. However, for the four sources with "intermediate X-ray luminosities" 
(B2 0755+37, 3C 31, 3C 317, and B2 0055+30; Lx = (2.4-5.2) x IQ-^LEdd), their X-ray origin 
is complicated. The X-ray spectrum of B2 0755+37 is dominated by the jet, while for 3C 31 and 
B2 0055+30 spectra are dominated by the ADAF. For the other one (3C 317), the contributions 
of the ADAF and jet are comparable. For the three least luminous sources (3C 66B, 3C 449, 
and 3C 272.1) which have Lx = (6.8 x 10-^ - 1 x W-^)LEdd, their X -ray spectra are dominated 
by the jet. The X-ray emission of 3C 449 is also interpreted by the sum of a jet and an ADAF, 
which requires higher quality data to further constrain it. Our results are roughly consistent 
with the predictions of Yuan and Cui (2005) . The "intermediate luminosity" here in our sample 
corresponds to the critical luminosity in Yuan and Cui (2005). However, the former is about 10 
times higher than the latter. The value of the critical luminosity depends on the ratio of the 
mass loss rate in the jet to the mass accretion rate in the ADAF. This ratio is adopted in Yuan 
and Cui (2005) from fitting the data of a black hole X-ray binary— XTEJl 118+480. The current 



result indicates that the ratio in FR Is is about 10 times higher than in XTE J1118+480. One 
possible reason could be that systematically the black holes in FR Is are spinning more rapidly. 
We note that the critical luminosity (or Eddington ratio) in this work is for the transition of 
X-ray emission dominated by the ADAF or by the jet, which is different from the critical value 
for the transition of the jet power dominated systems and accretion power dominated systems 
(Wu and Cao 2008 for more details and references therein). 

The kinetic luminosity, Lkin = — l)Afjetc^) can be derived from our modeling results. 

We use ryjet = L]^in/ M{\^Rs)c^ to describe the efficiency of the jet power converted from the 
accretion power, where M{lQRs) is mass accretion rate at 10i?5. We find that r^jet = 0.03 — 0.44 
for the sources in this sample (see Table 1), and most of them (six of eight) have r^jet significantly 
higher than 0.057, which is the largest available accretion energy at the innermost stable circular 
orbit (ISCO) for a non-rotating black hole. It implies that the black holes in these sources are 
spinning rapidly (so ISCO is smaller thus more accretion energy is available). 
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